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Abstract 

We present here a new MC study of ISB at finite temperature in a 
Z2 x Z2 A(/> 4 model in four dimensions. The results of our simulations, 
even if not conclusive, are favourable to ISB. Detection of the effect 
required measuring some critical couplings with six-digits precision, a 
level of accuracy that could be achieved only by a careful use of FSS 
techniques. The gap equations for the Debye masses, resulting from 
the resummation of the ring diagrams, seem to provide a qualitatively 
correct description of the data, while the simple one-loop formulae 
appear to be inadequate. 
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1 Introduction 



In 1974, in his classic paper on thermal field theory, Steven Weinberg, quot- 
ing S. Coleman, reported on the unusual behavior of certain particle theory 
models at high temperatures. Using the example of a simple 0(N\) x O(A^) 
scalar 4 model, he showed that it is possible to induce spontaneous sym- 
metry breaking at arbitrarily high temperatures. He named this funny phe- 
nomenon Inverse Symmetry Breaking (ISB) or Symmetry-Non-Restoration 
(SNR), depending on whether the ground state was disordered or ordered 
at zero temperature. In order to make this seemingly paradoxical state- 
ment easier to accept, he observed that in Nature there is a substance, the 
Rochelle salt, that does exhibit this sort of inverse behavior. Its crystals are 
orthorhombic below the lower Curie point (-18 °C), and monoclinic above 
it, which means that the crystal has a larger symmetry group at the lower 
temperature! Of course, Weinberg's statement was much stronger than this 
example, for he claimed that order could persist for ever, no matter how high 
the temperature, while the Rochelle salt, if heated enough, eventually melts. 

In any case, ISB and SNR would have probably remained relegated among 
the amenities of field theory, if it were not for the fact that this mechanism 
can be implemented in realistic particle models and then applied to important 
cosmological issues related to the thermal history of the Early Universe. This 
is not the right place to discuss these very interesting applications of ISB 
and SNR, but to give the reader an idea of the number and relevance of the 
issues that have been treated , we just make a list and quote some references. 
Following the historical order, ISB and SNR were applied to: the breaking of 
CP symmetry p|, the monopole problem ||, baryogenesis |§], inflation [§], 
the breaking of P, strong CP and Peccei-Quinn symmetries || and finally 
super symmetry 0. 

The potential importance of ISB and SNR for cosmology induced sev- 
eral authors to study these phenomena more accurately, than was done in 
Weinberg's original paper. His analysis relied, in fact, on a simple one- loop 
computation, and so it was natural to question if higher order corrections 
would have changed the scenario. This was not an idle question for at least 
two reasons. On one side, it is well known that at high temperature the 
reliability of perturbation theory is questionable, because powers of the tem- 
perature can compensate for powers of the coupling constants, spoiling the 
expansion. On the other side, and maybe more importantly, when one tries 
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to induce ISB or SNR in realistic particle models, one is generally forced to 
consider rather large couplings in the scalar sector, in order to overcome the 
action of the gauge fields, which always work in favour of symmetry resto- 
ration at high temperature, and thus it becomes important to explore the 
robustness of ISB and SNR with respect to the coupling strengths. 

The problem was attacked using a variety of approaches: large N expan- 
sions ||, the Gaussian effective potential ||, the constraint effective potential 
1(J, gap equations [pi , the average effective action |L2] and recently chiral 



effective Lagrangians [ 13j . The results were contradictory: some authors con- 
cluded that ISB and SNR are artifacts of the one-loop approximation, while 
others confirmed the existence of these phenomena, but observed that higher 
order-corrections have the effect of reducing the size of the parameter region 
for which they can occur. 

This state of things motivated us to start a program of MC simulations to 
study ISB, which seemed to us the only way to carry a fully non-perturbative 
analysis of this problem. We chose the simplest model that can exhibit ISB, 
namely a Z 2 x Z 2 two-scalars <fi\ theory. In the first run, we simulated the 



model in three dimensions |L4J and afterwards we followed with four ||15|| . In 
both cases we found no sign of ISB, despite the fact that, for the values of 
the coupling constants that were simulated, the one-loop conditions for ISB 
were strongly satisfied. Recently, the authors of ref. |]16| examined a^x 0(4) 
two-scalar 4 model, in four dimensions, and claimed to have found a clear 
evidence that ISB was taking place. 

In this paper we present the results of a new series of simulations of our 
Z 2 x Z 2 model. The strategy that we have followed is essentially the same of 



our previous paper flip] . Our two scalar model depends on five parameters, 
three coupling constants Ai,A2 and A, and two hopping parameters Ki and 
K 2 . We fixed once and for all the values of the coupling constants: with 



respect to [15], this time we enforced the perturbative conditions for ISB (in 



the direction of the field <p2) much more strongly and at the same time we 
took smaller values in order to be closer to the perturbative region. We then 
studied the phase diagram as a function of the two hopping parameters. In 
the {k\,K2) plane we found transition lines of first and second order, but 
since we were interested in ISB, we focused our efforts on the critical line for 
the breaking of the field 02- This line is roughly parallel to the Ki axis and 
the aim was to determine in what direction it shifts, when the temperature 
is increased: ISB requires that, for T > 0, it shifts towards smaller values of 



3 



K2, at least in a neighbourhood of the scaling region. 

As it is well known, finite temperatures are simulated by lattices with 
a finite extension N t in one direction, the temperature being T = l/(N t a), 
with a the lattice spacing. What made the simulations very hard is the 
extreme smallness of the effect that we show in this paper: even at the 
highest temperature, which corresponds to N t = 2, in order to detect the 
shift reliably, it was necessary to measure the critical values of K2, both for 
T = and for T > 0, with six significant digits! As it is well known, the only 
safe way to achieve such a huge precision is via an accurate analysis of Finite 
Size Scaling (FSS). This required large lattices and tremendous statistics: in 
one case, for example, we simulated a 20 4 lattice, with 8 x 10 6 iterations. 

For Nt = 2, the direction of the shift was favourable to ISB, but we cannot 
consider this result as conclusive. Due to the difficulty of the measurements, 
we could not simulate larger values of Nt, as it is necessary in order to make 
sure, via scaling analysis, that things would go in the same way for all A^'s. 

As a check, we compared the theoretical predictions with the MC value 
of the critical temperature for which we observed ISB. While the simple one- 
loop estimate is grossly incorrect, we found a reasonable agreement with the 
value obtained from the gap equations, which result from resumming the 



ring diagrams of the perturbative series [pT |. These equations seem to give 



a qualitatively correct description of the MC data, and explain as well why 
we did not find ISB in our previous simulations in |15| . 

We close this introductory Section with a plan of the paper: in Section 
2 we introduce our lattice model and discuss its phase diagram at T = 0. 
In Section 3 we review the perturbative picture of ISB, while in Section 4 
we discuss in detail our strategy to study ISB on the lattice. In Section 5 
we discuss how the results of the simulations compare with the theoretical 
predictions, while Section 6 contains a discussion of the recent MC studies 
15H and [16]. Finally, Section 7 contains our concluding remarks. 



m 
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2 The model and its lattice formulation 



We consider the theory for two real scalar fields in 4 euclidean dimensions, 
described by the bare (euclidean) action: 

s=f**{T, [~(^#V~™f )2 <^ (i) 



i=l,2 

What will be essential for ISB and SNR, in the above action the quartic bare 
coupling can be negative. If g^ > the condition of boundedness from 
below of the potential is satisfied if: 

si 0) >o, #>0, 

9? > 9 g^. (2) 

When regularized on an infinite four- dimensional cubic lattice of points Q 
with lattice spacing a, the above action is replaced by its discretized version 



(o) 

r(°)2,f,2 f \ , 9i ^4 



+ 



+ 9 — (3) 

where AJj 1 ' is the lattice derivative operator in the direction fx: 

AW^(x) = <M* + /*Q-<M*) . (4) 

We find it convenient to measure all dimensionful quantities in (H) in units 
of the lattice spacing; thus we define: 

4>i,L( x ) = a ®iA X )' m i = a2m t )2 ■ (5) 

In terms of the dimensionless quantities the lattice action now reads: 

1 1 ( 0) 

Sl = E { E [Ei( A ?W) 2 (x) + ^rnf )2 4>l L {x) + ^-< L (ar) 
xen i=i,2 a* ~ 



+ 
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The standard lattice notation is obtained with a further redefinition of the 
fields and couplings in (H) according to: 

(piA x ) = VKi<f>i, r , gf ] = ™f )2 = 2- — — — — (7) 

and 

9 W = — - (8) 
After these redefinitions we get our final form of the lattice action 



Sl = { [ ~ K i <f>i,r<f>i,r+fl + M0i,r ~ + 



> 2 



reZ 4 i=l,2 



+ 



+ WlAlr}- (9) 



For generic values of the parameters, the action has a Z2 x Z2 symmetry 
which can be spontaneously broken. The model is expected to have only one 
fixed point, the Gaussian one, corresponding to mj[ = = 9i = 9% = 
g(°) = (or K\ = K2 = 1/4 and Ai = A2 = A = 0), which has a null attractive 
domain in the infrared. In 4 — e dimensions there exist five more fixed points: 
1) the Heisenberg fixed point, for a e g^ = a e g^ = 3a e g^ = g* H and m^ 2 = 



m 



(0)2 



m 2 ^ < 0, where the symmetry of the model is enhanced to 0(2). 



2) the three Ising fixed points, for g^> = and a € gl°^ = gf, mf^ 2 = m|* for 
some i, where the model splits into two independent 4>\ models. 

3) the Cubic fixed point, for a e g^ = a e g^ = gJ/2, mf^ 2 = raj*, which again 
splits into two independent <p\ models, after a 7r/4-rotation of the fields. 

The phase diagram of the model @ (at T = 0), for fixed values of Ai, A 2 
and A, is shown in Fig. 1. 

There are four distinct phases separated by four critical lines A, B, D and 
E: the disordered one corresponds to the lower left corner of the picture, while 
the upper right corner represents the totally ordered phase. The wedges in 
between them are partially ordered phases: in the left wedge the field 02 is 
ordered, while <f>\ is disordered; in the right wedge it occurs the contrary. 

If there was no interaction between 0i and 02, we would have two inde- 
pendent 4 models and then the critical lines would all be of second order. 
The lines A and E would form a straight line parallel to the K\ axis, and 
similarly the lines B and D would form a straight line parallel to the K2 axis. 
The presence of the interaction bends them as can be seen in the Figure. We 
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Figure 1: Phase diagram at Ai = 0.3375, A 2 = 0.01125, A = -0.112. 



managed to determine the order of the critical lines only at some distance 
from the region C, where the four lines merge together. The critical line A is 
the only one that we studied accurately: away from C, it is of second order 
and we observed for it Gaussian critical exponents, as it had to be expected 
from triviality. As for the remaining critical lines, there is a clear evidence 
that D is of second order, while B and E appear to be of first order. Close 
to C, things become unclear: while on small lattices it becomes impossible 
to distinguish the various phase transitions, on large ones the onset of strong 
metastabilities prevented us from getting any indications at all. 

Understanding the phase diagram in the neighbourhood of C is of course 
extremely important. An important issue to address is whether the four 
critical lines actually meet at a single point, a double critical point that we 
will call C, and what is its order. 

Indeed, the indications coming from perturbation theory suggest that, 
when g R is negative, as required for ISB or SNR to occur, C, if it exists, 
should be a first order critical point. This can be seen by integrating the 
one-loop renormalization group equations (RGE), in the scaling limit of large 
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correlation lengths & for both fields 



& = q e , f > i, 

where constants. In this limit the RGE of our model read: 



dt 167T 2 

dg R 3 



dt 167T 2 1 

+ W, (10) 



dt 167T 2 

where t = log£. Numerical integration of these equations shows that, if 
one starts from a set of initial values for the renormalized couplings with 
g R negative, g R and g R are driven by g R towards negative values in a finite 
interval of t, something that is generally believed to indicate a first-order 
phase transition. The fact that for negative values for the coupling g R the 



double critical point really is of first order, is shown in ref.[fL8fl, where the 



phase diagram of our model is studied, using a suitable truncation of the exact 
non-perturbative flow equations, for the particular case gi = gf \ mf' 2 = 
m^ 2 1 when an extra 0i «-> 02 symmetry is present (this case is not directly 
related to ISB or SNR because this extra symmetry excludes a priori the 
possibility of these phenomena). 

If this picture of the double phase transition for g R < turns out to be 
correct, one may question whether in this case such a model can be used at 
all as an elementary particles theory. The answer to this question depends 
on the actual strength of the first-order phase transition. The fact that the 
correlation lengths of one or both fields remain finite near C is not necessarily 
a problem per se. In fact, even if C were of second order for both fields, 
triviality would prevent one from taking the correlation lengths to infinity 
anyway. Thus, the relevant question is to see if the correlation lengths become 
sufficiently large near C, for the presence of the lattice cut-off to become 
negligible in low-energy amplitudes. This in turn depends on the absolute 
magnitude of the coupling constants. For small couplings, C should be only 
weakly first order and one should be able to achieve large correlation lengths 
for both fields. In realistic particle models for ISB and SNR the problem 
may be more serious: the presence of gauge couplings, that work in favour 
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of symmetry restoration, usually forces one to take large negative couplings 
g R , for ISB or SNR to occur |[L1]| , and this can then give rise to rather strong 
first order phase transitions in the scalar sector, thus spoiling the model. 



3 High-T perturbation theory 



In this Section we briefly review the predictions of perturbation theory on the 
high-temperature behavior of the model ([!]). We shall see that, for certain 
choices of the couplings, the symmetric vacuum of the theory becomes unsta- 
ble at sufficiently high temperatures. In order to explore this phenomenon, 
we need to compute the leading high-T contribution to the second deriva- 
tives of the effective potential in the origin, or equivalently the 1PI 2-point 
functions for zero external momenta p. A standard one-loop computation (in 
the imaginary time formalism) |19| in our model gives the result: 



M'f(T) = -T ( ?\p = 0,T) = mi + 



M|(T) 



-lf ) (p = 0,T)=m^ + ^ 



rp2 

24 
T 2 



T 



T J 



9 ?ft(f) + <A(^ 



,(n) 



where mf represent the renormalized masses at T = and h{x) is the func- 
tion: 



6 f°° 

h ( x ) = ~2 / d y 

7T Z JO 



It is clear that h(x) is positive, monotonically decreasing and that it ap- 
proaches zero when x — > oo. For small values of x, h(x) has the asymptotic 
expansion: 

h(x) = 1 - -x - -^x 2 ( In ^- + 7 - \ V (13) 



7T 



7 being the Euler constant. In the limit of very high-temperatures rrii/T 1, 
eqs.fllll) reduce to: 



M 2 (T)=m 2 2 + ^(g? + g 



(14) 
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Assuming for simplicity that the bare couplings are so small that they can 
be identified with the renormalized ones, and recalling that the coupling 
g(o) gR can kg negative, it is easy to see that for g R /g R sufficiently large, 
the range (§) of stability for the potential includes values of g R such that, 
say, g R + g R is negative (it can be proven easily form eqs.(Q) that one cannot 
have g R + g R < at the same time). In this case we see from eq.(|14]) that 
M|(T) becomes negative at sufficiently high temperatures, irrespective of its 
value at T = 0. This is the essence of the phenomena of SNR and ISB: in 
multiscalar models some Debye masses can become negative and so one can 
have spontaneous symmetry breaking at arbitrarily high temperatures. We 
now focuse on ISB, the phenomenon to be explored in this paper: in this 
case we take m\ > 0, m\ > 0, so that the vacuum is disordered at T = 0, 
and again assume that g R + g R < 0. If one starts increasing the temperature, 
it will happen that, at a certain critical temperature T c , the field 02 will 
undergo spontaneous symmetry breaking and will stay broken at all higher 
temperatures. The field 0i, instead, will remain disordered. The critical 
temperature T c can be estimated from eq. (|T4|) to be: 



T c I 24 



™2 V 9 + 92 



(15) 



Even though they capture the main features of ISB and SNR, at a closer 
inspection, eqs.(|14D and ( |15"D appear unsatisfactory in two respects. First of 
all, the condition for the high-T instability in, say, the 02 direction, g R +g R < 
0, as deduced from the second of eqs.([L4]), does not involve g R . Since the 
bound (0), together with the condition g R + g R < 0, typically implies a 
hierarchy among the couplings g R ^> \g R \ > g R , it is clear that next-to- 
leading order corrections involving g R can produce significant corrections to 
the lowest order result. The second fault of eqs.(|l4]) and (|l"5"D is the most 
important for our study of ISB. It is the fact that the estimate ([15]) of the 
critical temperature does not involve m\: this looks physically incorrect, for 
one can imagine that for m\ ^> vri2 the field 0i should decouple from 02 for 
temperatures T<mi and so ISB or SNR should not occur before the scale 
mi is reached. Both limitations of the one-loop picture can be overcome by 
performing the resummation of the so called ring diagrams. 

As it is well known [I7|, these diagrams represent the dominant next-to- 



leading corrections at high temperatures and their inclusion leads to a more 
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accurate description of both the phase transition and the high temperature 
behavior of the thermal masses. This infinite resummation leads to self- 
consistent gap equations for the Debye masses Mf(T), which in our case 
read 



.1, ■-} 



|t + ff M*2) + f^(*i) , (16) 



where we defined X, = Mj(T)/T. For any given temperature T, the symmet- 
ric vacuum is stable if the above equations admit real (positive) solutions 
for X{. At sufficiently low temperatures, this is clearly the case if we take 
m\ > 0, m\ > 0. One can now gradually increase T and follow the evolu- 
tion of Xi and x<i'- if one of them, say x^, approaches zero at T = T c , and 
if for T > T c the gap equations do not admit anymore real solutions, one 
can argue that for T > T c there will be ISB in the </>2 direction. It is clear 
that this can happen only if g is negative and its absolute value is large 
enough in comparison with g%. A detailed study of the conditions required 



is given in the first of refs.[|ll[], where it is shown that the parameter region 
for which ISB or SNR occur is indeed smaller than that predicted by the 
naive eqs.(jTJ|). An analogous result was later confirmed by studies based on 



different non-perturbative techniques [12 



Assuming now that at high temperatures there develops an instability in 
the 02 direction, a better estimate of the critical temperature T c , than that 
given by the one-loop equations (|T5D, can be obtained by setting x 2 = in 
eqs. (|16D and then solving the system: 

A = ^| + y a {9iK^) + g R ) , 

= ^ + ^ + , (17) 



where we have used h(0) = 1. Upon comparing the second of eqs. (|T7|) with 
the second of eqs. (|l"4|) , it is clear that the inclusion of the ring diagrams 
pushes the phase transition to higher temperatures. This is so because the 
negative coupling g , which is the cause of ISB, comes with the factor h(x\), 
which is always less than one. Since X\ is larger then mf/T c 2 , as it is clear 
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from the first of eqs.flTTD, we see that a large mass m\ will make ISB harder 
to achieve. In fact, the influence of m\ can be so strong as to cause the 
disappearance of the symmetry-breaking phase transition predicted by the 
one- loop formula. As we shall explain in greater detail in Sec. 6, this effect 



probably explains why we did not observe ISB in our previous work [15 



4 ISB on the lattice 

We saw in the previous Section that, for + g R sufficiently negative, con- 
tinuum perturbation theory predicts ISB for the <p2 field. Even if the mass 
parameters are chosen in such a way that the ground state is disordered at 
T = 0, the field 02 should develop a non- vanishing vev above some critical 
temperature T c and should remain ordered at all higher temperatures. The 
field 0i, instead, should remain disordered. 

In principle it is very simple to study ISB on the lattice. As it is well 
known, finite temperatures are simulated by lattices with a finite extension 
in the "time" direction (and infinite extension in the remaining "space" di- 
rections, as required by the thermodynamic limit), the temperature being 
related to the number N t of sites by the relation: 

(18) 



N,a 



with a the lattice spacing. So, in order to study ISB, one should simply 
check how the critical lines of the T = phase diagram, shown in Fig.l, 
shift in the (k\, k^) plane (for fixed values of Ai, A2 and A) as a function of 
the temperature, which means concretely as a function of N t (in principle, 
one could vary the temperature also by varying the lattice spacing in the 
time direction, but we have not used this method). Since we are specifically 
interested in the behavior of the field 2 , it is sufficient to focuse our attention 
on the critical line of this field, the line A-E of Fig.l. If, for the selected values 
of Ai, A 2 and A, ISB occurs, it should happen that, when the temperature 
is increased, namely when N± is diminished, the 02-critical line penetrates 
more and more deeply in the disordered region of the T = phase diagram, 
at least in a neighbourhood of the double critical point C. The effect should 
be clearest at the highest temperatures, i.e. for the smallest values of N t . To 
help the reader visualize the situation, we have shown in Fig. 2 the expected 
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Figure 2: Expected (qualitative) temperature-behavior of the critical line for 
the 02 transition, in case of ISB. The region studied by us is the one to the 
left of C. 

iV r behavior of the 2 -critical line, in case of ISB. Some comments are in 
order. 

We expect that far from the double critical point C the N t behavior of 
the 02-critical line should always be the normal one, namely the critical 
line should move towards greater values of k 2 when N t is diminished. This 
conjecture looks reasonable, if one considers that far from C, the field 0i 
should be far from critical and its correlation length should become very 
small. This implies that the influence of the field 0i over the field 02, near 
the critical line of the latter, should be negligible and so the behavior of 02 
should be that of a normal 4 theory, where an increase of temperature shifts 
the transition upwards. This means that ISB should be visible, if it happens 
at all, only close enough to C, in the scaling region where the correlation 
lengths of both fields become large. This explains why, in Fig. 2, we have 
drawn the critical lines for N t = 2,3,4,5 below that of the T = theory 
(Nt = oo) only near C, while away from C they all lie above it. 
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Such a behavior poses a very serious problem: imagine that for the few 
values of N t simulated, one has found a result like that shown in Fig. 2, 
which supports ISB. How can one convince oneself that nothing will change 
for greater values of N t l One could suspect that for N t large enough the 
critical line might pass entirely above that of the T = theory, thus leading 
to the disappearance of ISB. In practice, it is then very important to have a 
criterion to decide if, whatever behavior is observed with the few iV t 's that 
are simulated, it will not change when going to higher iV t 's. The common 
practice to address this question is to search for scaling in the N t dependence 



of some observable: this is what we did in [|15[], where, for a fixed value of 
K\ and several values of Nt, we measured the critical values of k 2 , K^iN). 
The values we got were all greater than the critical value k% for T = 0, a 
result against ISB (very likely, the value of K\ used for the simulations in 
[ |T5| was too far from C, for ISB to be visible. We were probably in the 
region of Fig. 2 on the left of C, where the behavior of the 02 critical line 
becomes the "normal" one. We shall have more to say on this in Sec. 6.) 
. The good scaling of k^N) with Nt convinced us that no change would 
have occurred for larger N t . We think that applying this method to make 
sure that ISB survives for all values of N t would be very difficult, because 
we believe that in order to observe ISB for a number of values of N t large 
enough for a scaling analysis to be possible, one needs look very close to the 
critical point C, where the large correlation lengths make it necessary to use 
very large lattices for a reliable measurement of the critical couplings. 

As a matter of fact, in the new series of simulations presented here, we 
have been able to simulate only one value of N t , N t = 2, and we needed a 
very large statistics in order to detect reliably the shift of the 02 critical line. 

Let us briefly explain the method we have followed, postponing the details 
to the next two subsections. 

The first part of the method is exactly the same as that of ref.[[lfj, but 
it is useful to review it here again. We fixed once and for all a set of values 
for A 1; A 2 and A. This time, we took A x = 0.16875, A 2 = 0.0001125, A = 



—0.00784. There are two important differences with ref. [15]: 

a) the new values of the couplings are smaller than the old ones |15[] . This was 
done to be closer to the perturbative region, but it increased tremendously 
the difficulty of the simulations, because the shift of the critical line with N t 
was this time much smaller than before; 

b) more importantly, the ratio | A| /A2 is now much greater than before (more 
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N t 


N s 


# runs 


Thermalization sweeps 


# sweeps 


2 


6 


1 


5 x 10 4 


4 x 10 6 


2 


8 


1 


10 5 


8 x 10 6 


2 


10 


1 


2 x 10 5 


10 7 


2 


12 


1 


3 x 10 5 


10 7 


2 


16 


2 


2 x 10 5 


8 x 10 6 


2 


20 


4 


8 x 10 4 


6.4 x 10 6 


6 


6 


1 


5 x 10 4 


4 x 10 6 


8 


8 


1 


10 5 


8 x 10 6 


10 


10 


1 


2 x 10 5 


4 x 10 6 


12 


12 


1 


3 x 10 5 


8 x 10 6 


16 


16 


4 


2 x 10 5 


9.6 x 10 6 


20 


20 


6 


2.4 x 10 5 


9.6 x 10 6 



Table 1: Statistics for the simulations on lattices of sizes N t x iV s 3 . 

or less 70 against 10), and so there is now a much stronger push towards ISB. 

After a quick analysis of the phase diagram at T = 0, we selected a 
value of Ki, Ki = 0.24, to the left of C, and then we measured accurately 
the corresponding critical value k 2 of k 2 along the (f> 2 critical line, always at 
T = 0. 

Near this point, we measured the T = values of the renormalized cou- 
plings and the correlation lengths of both fields, for the same value of Ki, but 
for a value of k 2 slightly less than k,^, in the disordered phase (the results of 
these measurements are given in Table 2). The value of £i ~ 1 that we got 
is not very large, but still acceptable for us to believe that we were probing 
the continuum region. In any case, we could not achieve larger values for 
£i, because near to C we encountered strong metastabilities that made any 
simulations impossible. We notice also that the value of g R is quite large 
and not very accurate; on the contrary, what we think is the essential, we 
see that g R and \g R \ are both very small and that the condition for ISB is 
implemented very strongly, since \g R \/g 2 — 10. The one-loop formulae obvi- 
ously predict that ISB should occur for these values of the couplings, and it 
can be checked that the ratio \g R \/g 2 is large enough for ISB to survive also 
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Asymmetric Lattice. Binder Cumulants 




Figure 3: Binder cumulants for the 2 x iV s 3 lattices (N s = 6, 8, 10, 12, 16, 20) 
at Ai = 0.16875, A 2 = 0.0001125, A = -0.00784 and m = 0.24. 

the inclusion of the ring diagrams. This convinced us that we had found a 
good simulation region. 

Having performed all these measurements at T = 0, we turned to finite 
values of N t . Due to the extreme smallness of the effect, we were able to detect 
the shift of k c 2 (for K\ = 0.24) only for the smallest value of N t , N t = 2. Even 
in this case, the shift was of the order of one part in hundred thousands and 
so larger values on N t look out of reach. Since this shift turned out to be 
towards smaller values of k 2) our result can be considered as favourable to 
ISB. 

In the next two Subsections, we describe the details of the numerical 
simulations and how we managed to measure the critical values of k 2 with 
the high accuracy that was needed. 
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4.1 MC Simulation: determination of k\ 

For this new simulation, we used the Metropolis Algorithm (cluster was not 
efficient), with 3 hits and acceptance around 60%. We simulated lattices 
of sizes 2 x and Nf, with N s = 6, 8, 10, 12, 16, 20. In the new region of 
parameters the field fluctuations are strong and the transitions for different 
lattice sizes are very close. We had to accumulate very large statistics in 
order to obtain accurate results. In Table p] we summarize the statistics 
for every lattice: number of runs, sweeps left for thermalization and total 
number of sweeps for every run. 

As we said above, all our simulations were performed for a single value 
of Ki, Ri=0.24. A fast simulation on the 6 4 lattice gave us the approximate 
position, k 2 = 0.24903, of the phase transition of the field 02 (the field 
0! remains disordered across the transition). All the subsequent massive 
simulations were performed at this point, and the results were extrapolated 
in a narrow ^-interval around it by means of the Spectral Density Method 
[p2]l (SDM). The observable used to locate with high precision the phase 
transition was the Binder cumulant: 

U i (N.,K 1 ,K 2 ) = --±^, (19) 

where stands for the magnetization of the field 0j. 

FigsJ^ and ^ show the values of U2 obtained by extrapolating the results 
of the simulations by means of the SDM. Every pair of curves (N s and N s ) 
determines a crossing point, that gives an estimate K^Ng, N s ) of the critical 
point k%. These values do not suffer from finite size effects, except for correc- 
tions to scaling, which can be parametrized and used to obtain k^, according 
to the scaling law: 

1 — 

k* 2 (N s , bN s ) - 4 = W —Nr^ lfv , (20) 

where u is the exponent for the corrections to scaling. In the calculation of 
the error, due account was taken of the fact that the pair crossings were not 
all independent of each other (out of the iV(iV — l)/2 crossings between N 
curves, only 2N — 3 are independent of each other). 

As for the exponents occurring in eq. (P0|), we had to distinguish the 
symmetric {Nf) lattices from the asymmetric (2 x N^) ones. In the case of 
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Symmetric Lattice. Binder Cumulants 




0.24900 0.24901 0.24902 0.24903 0.24904 

Figure 4: Binder cumulants for the Nf lattices (N s = 6, 8, 10, 12, 16, 20) at 
Ai = 0.16875, A 2 = 0.0001125, A = -0.00784 and Ki = 0.24. 

the latter, since the scaling parameter is N s , while N t is fixed, according to the 
hypothesis of dimensional reduction and universality, we used the exponents 
of the Ising model in three dimensions, namely v = 0.63 and uj = 0.8. As 
a check, we also computed v directly from the data obtaining values fully 
compatible with the above one (see Fig. [7]). 

In the case of the symmetric lattices, we used the mean field exponents 
uj = 0, v = 1/2. As seen from eq.(pO|), these values imply that, apart from 
logarithmic corrections, all the crossings should occur for the same value of 
K2, and this is approximately what we found. A better estimate of k 2 was 
then obtained by setting u> ^ in eq. (|20|) and taking the limit, for uj — > 0, 
of the corresponding values of k\. In the process, we observed a very fast 
convergence. 

The result of the extrapolation for N s — > oo is shown in Fig. |5|. We got 
K c 2 (N t = 2) = 0.249005(2) and k c 2 (T = 0) = 0.249022(2). This means that 
the disordered region diminishes when T is increased, as ISB requires, at least 
for the case of the highest temperature reachable on the lattice (Nt = 2). 
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0.24904 



0.24902 



0.24900 — 




0.24898 — 



0.24896 



0.000 0.001 0.002 0.003 0.004 0.005 

N s -"- lA, (i-b-")/(b 1/y -i) 



Figure 5: Fits to obtain k 2 from the crossings of the Binder cumulants. The 
results are re*(T = 0) = 0.249022(2) (filled) and K c 2 (N t = 2) = 0.249005(2) 
(empty). 



4.2 MC Simulation: maximum of the Binder cumulant 
derivative 

Due to the extreme smallness of the shift in the value of k 2 , we found it 
opportune to compute the two values of k 2 also in another way, completely 
independent on the previous one. Indeed, for every value of N s , one can 
determine the point k 2 *(N s ) where the derivative dUijdKi reaches its maxi- 
mum value, namely C v . It is well known that, in the limit N s — > oo, k 2 *(N s ) 
approaches k 2 according to the scaling law (including corrections to scaling): 

k* 2 *(N s ) - k% oc JV^i _ an;") (21) 

which allows us to compute k 2 . Fig.|6] shows the corresponding fits. In this 
computation, we used the same critical exponents discussed in the previous 
subsection. We obtained K c 2 (N t = 2) = 0.248997(11) and k c 2 (T = 0) = 
0.249017(2), which are fully compatible with the previous results. 
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0.24910 



0.24908 



O 

0) 



U 0.24906 



0.24904 



0.24902 



7.5 



• N t =N s /e(°°) = 0. 249017(2) 
o N t = 2 /c(°°) = 0. 248997(11) 




20.0 



Figure 6: Fits to obtain from the maximum of the Binder cumulant deriva- 
tive. (The exponents used were v = 0.63, uj = 0.8 for N t = 2 and v = 1/2 for 
N t = N K ). 



As a final check, we computed v from the maximum of dUi/dn2, C v , 
which diverges with the lattice size as gain, we obtained consistent 

results, as can be seen in Fig. |7[ 

In conclusion, this time, differently from what we reported in [15 1, we 
found that, when N t = 2, the phase transition of the field 
value of «2 which is smaller than at T = 0. 
long time, because, for the small values of A2 and A used now, the shift of 
the critical </>2-lme is extremely small and we had to measure both values of 
k c 2 with six significant digits, a job that took us months of computer time! 

So, the results of our new simulations are favourable to ISB, but cannot 
be considered as conclusive because they refer to a single value of N t . 



occurs for a 
The two measures of k c 2 took a 



20 



30000 




Figure 7: Fits to obtain v by using the maximum of the Binder cumulant 
derivative. 



5 Comparison with Perturbation theory 

As we said in Section 3, simulating only a finite number of N t values is 
per se not sufficient to claim ISB, for one cannot rule out the possibility 
that the phenomenon would disappear for larger values of N t . In fact, the 
case here is even worse, since we have at our disposal just one value of 
Nt, Nt — 2. In order to increase confidence in our result, we attempted a 
comparison of our MC data with the perturbative predictions. While the 



one-loop formulae, eqs.(|T^), lead to an estimate of the critical temperature 
far from the MC value, we found a reasonable agreement with the predictions 
of the gap equations, eqs.([T7|). 

We shall now explain how this comparison was done. Let p = (ki,/^) 
be the point of the 02-critical line for N t = 2, that we identified in our 
simulations. As we saw, p belongs to the disordered phase of the T = 
phase diagram. If we knew the T = values of the renormalized couplings 
and correlation lengths £i and £2 (we measure them in units of the lattice 
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N s 


6 


6 




02* 




6 


0.9581(2) 


10.78(1) 


26(6) 


0.0261(9) 


-0.35(1) 


8 


0.9528(3) 


11.48(2) 


41(11) 


0.0296(15) 


-0.29(1) 



Table 2: T = values of the correlation lengths and renormalized couplings 
for ki = «i = 0.24, k 2 = k 2 = 0.2488 (lattice sizes N*). 



spacing) at p, we could get an estimate of the critical inverse temperature 



N£ by solving eq.(|T3) or eqs.(|T 

If perturbation theory works, we should get for N% a value close to two. 

This procedure, even if simple from the conceptual point of view, is diffi- 
cult to carry out: since p is very close to the critical line of the T = phase 
diagram, the correlation length £2 is huge there (in our case, we estimate it 
to be larger than 40). Thus a measure of the T = renormalized couplings 
and correlation lengths at p requires 4d-lattices with iV s > 40, which is just 
too much for computers available today. 

In order to avoid the difficulties connected with exceedingly large correla- 
tion lengths for the field <j) 2 , we measured the renormalized couplings and the 
correlation lengths (for T = 0) at a point p = R 2 ), having the same K\ 
coordinate as p, but with R 2 slightly less than k 2 , an d such that £ 2 was more 
or less ten (see Table |2|). Because of time limitations, we did not carry out the 
thermodynamic limit in this series of measures and contented ourselves with 
simple estimates on rather small lattices. The quantitative considerations 
that follow have to be taken, then, as simply indicative. 

Afterwards, we computed the values of #f , g^, gn and of £1 and £ 2 at 
p, by integrating a set of (T = 0) RGE's in the interval R 2 < k 2 < k 2 . We 
evaluated the /3-functions to two-loops order, following a procedure strictly 
analogous to that discussed in p0[ . Of course in our case there was a larger 



number of (3 functions to be computed. Since the value of £1 at p was only 
roughly one, we included in our /5-functions the corrections to scaling upto 
one-loop order. 

The reader might question the reliability of these perturbative RGE's, 
since the value of gf at p is rather large and quite uncertain. We believe that 
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this is not a problem, for two reasons. The first one is that the integration 
interval in K2 is very small. The second, and most important, is less obvious: 
it turns out that g R does not appear in any of the (3 functions (at least to the 
two-loops order that we have examined). Its absence can be understood by 
considering that, since we are moving in the vertical direction in the (k\, K2) 
plane, £1 remains finite (and practically constant). 

After the numerical integration, it turned out that the running 
g R and £1 in the tiny K2 interval considered was completely negligible, while 
£2 varied from 11 to something like 45. Since the final value of £2 was very 
sensitive to its initial value, (and practically insensitive to the initial values 
of the other parameters), we took advantage of our accurate knowledge of 
k\{T = 0), and tuned the initial value of £2 in such a way that K2 would 
approach the MC value of k 2 (T = 0) in the limit £2 ~^ 00. The initial value 
of £2 that we obtained in this way, £ 2 — 11.85(5), is very close to the MC 
one, as can be seen from Table 2. The value of £2 a t P turned out to be 
&(p) — 45(5): the error is basically due to the uncertainty in the MC values 
of k c 2 (T = 0) and n c 2 (N t = 2). 

We thus solved eq. (p~5|) and eqs. (|i"7|) with respect to N£, using as inputs 
the values of gf , g R , g R and £1 given in Table 2 and the above value of &(p)- 
In this way we obtained the following estimates of iV t c : 

^(one-loop) = 4g(6) ? ^(gap) = ^ 

It is clear that ]\/" t c ( ono ~ loo p) j s f ar f rom MC value N t = 2, while N^ 9 " 1 ^ is 
quite close. It seems, then, that resummed perturbation theory provides a 
rather accurate description of our data. 



6 Some comments about recent results on 
ISB on the Lattice 

Recently, there have been two other MC studies of ISB in four dimensions, 



one by our group |To] and another by Jansen and Laine [JTHJ] . While in the first 
one no sign of ISB was found, the latter authors claimed to have seen clear 
evidence of ISB and good agreement with resummed perturbation theory. 
We would like to make here a few comments on these works. 
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We shall start from our own previous work. The method followed there 
was very similar to that used in this paper, but the simulations were per- 
formed for a different choice of parameters. We think now that the essential 
reason for the negative result is the fact that the conditions for ISB were 
not enforced strongly enough. Near the simulation point, the absolute value 
of g R was roughly the double of the value of g R : according to the one loop 
formulae flUD , this was more than enough to induce ISB, but nevertheless 
we did not find it. The reason of the failure is hidden in the gap equations 
flloT) . Near our simulation point the correlation length of <f>\ was approxi- 
mately one. This implies that, even at highest temperature (corresponding 
to N t = 2) and so in the conditions most favourable to ISB, mi/T = Nt/£i 
was approximately equal to two. By looking at the gap equations, it is easy 
to convince oneself that, with a value of m\/T ~ 2, a \g R \ which is twice 
g R is not enough to induce ISB. The reason is the suppressing factor h{x\) 
in front of g R in the second of eqs.(|16|): since X\ m 2, and h(2) ~ 0.17 it is 
clear that one would have needed a \g R \ at least six times bigger than g R . In 
order to be on the safe side, in the new simulations presented in this paper, 
since £i was still of order one, we selected \g R \ ~ 10 g R which now was large 
enough. 

As for ref.[|l6j, the authors searched for ISB in an 0(4) x Z2 two scalar 
model. Their approach is very different from ours, and it is useful to briefly 
review it here. Using perturbation theory (with two-loops accuracy), they 
computed the relations giving the (bare) lattice parameters, in terms of the 
renormalized masses and couplings, in the symmetric phase of their model. 
Then, they performed a set of simulations on lattices with N t — 2, 4 (the 
few short runs on symmetric lattices, according to the authors, should not 
be given much importance) for a number of values of the bare couplings, 
corresponding to decreasing values for the renormalized mass mi of the field 
expected to undergo ISB (0i), and constant values for the remaining renor- 
malized parameters. In this way they were able to vary the ratio T/mi in a 
large interval, even holding N t fixed, or varying it very little, by just varying 
the renormalized mass m\. In the simulations they measured the expectation 
value of the modulus of 0i, and the results turned out to be fully compatible 
with the predictions of (resummed) perturbation theory. In particular, the 
field 0i appeared to be disordered for values of T/rrii small enough, and or- 
dered for larger values of that ratio, which is a sign of ISB. Compared with 
our approach, this procedure presents the advantage of avoiding the very de- 
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manding task of an accurate determination of the critical point, but on the 
other side it is necessarily limited to the region of very small couplings, in 
order for perturbation theory to be applicable. Indeed, in ref.|TjJ, perturba- 
tion theory was used not only to compute the relations between the bare and 
renormalised couplings, as we said above, but also to monitor the behavior 
of < |0i | > in the thermodynamical limit. In view of the long lasting debate, 
whether ISB and SNR are genuine effects or rather artifacts of perturbation 
theory, we preferred to pursue the approach presented in this paper which, 
even if much more expensive, in terms of computer time, does not make any 
recourse to perturbation theory. As a reward, we were able to probe much 
larger values of the coupling constants, something which is especially impor- 
tant in view of the applications of ISB to realistic cosmological scenarios, 
where large couplings in the scalar sector are usually needed in order to en- 
force ISB, as was pointed out at the end of Sec. 2. Moreover, we could also 
study in detail the nature of the symmetry-breaking high-T phase transition. 



7 Conclusions 

We have studied ISB at high temperature in four dimensions in a two-scalar 
4 model with Z 2 x Z 2 global symmetry. For fixed values of the coupling con- 
stants Ai,A 2 ,A and of the hopping parameter K\, we measured the critical 
value of the hopping parameter k 2 , for which the field 2 undergoes sponta- 
neous symmetry breaking. This measure was performed both at T = and 
at T > 0. It was found that the value of k\ at T > is smaller than that 
for T = 0, which is a signal of ISB. Due to the extreme smallness of the 
shift, it was necessary to measure the values of k\ with very high precision. 
This accuracy was reached via an accurate analysis of FSS on large lattices 
with very high statistics. The extreme difficulty of the measurements did 
not allow a scaling analysis of the dependence of the result on the number of 
sites N t in the time direction. For this reason, we cannot consider our results 
as conclusive. 

In order to increase our confidence in the result, we compared the MC 
value of the critical temperature T c with the theoretical value. While the 
simple one-loop estimate turned out to be grossly incorrect, we found a rea- 
sonable agreement with the value of T c obtained from the gap equations (|T7|) , 
which result from the resummation of the ring diagrams. 
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So, in conclusion, we have found some evidence that ISB is possible and 
that resummed perturbation theory provides a qualitatively correct descrip- 
tion of the phenomenon. Nevertheless, we think that our results are not 
conclusive and that further work is needed before ISB can be claimed with 
confidence. 
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